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CONFORMAL AND NONCONfORMAL surfaces 

Hydrodynamic lubrication is generally characterized by surfaces that 
are conformal. That is, the surfaces fit snugly into each other with a high 
degree of geometrical conformity, so that the load is carried over a rela- 
tively large area. Furthermore the load-carrying surface area remains 
essentially constant while the load is increased. Fluid-film journal and 
slider bearings are conformal surfaces. In journal bearings the radial 
clearance between the shaft and the bearing is typically one-thousandth of 
the shaft diameter; in slider bearings the inclination of the bearing sur- 
face to the runner is typically one part in a thousand. 

Many machine elements have contacting surfaces that do not conform to 
each other very well. The full burden of the load must then be carried by a 
very small contact area. In general the contact areas between nonconformal 
surfaces enlarge considerably with increasing load but are still small com- 
pared with the contact areas between conformal surfaces. Some examples of 
these nonconformal surfaces are mating gear teeth, cams and followers, and 
rolling-element bearings. 

The load per unit area in conformal bearings is relatively low, typi- 
cally only 1400 N/m^ and seldom over 3b 000 N/m^. By contrast, the load 
per unit area in nonconformal contacts, such as those that exist in ball 
bearings, will generally exceed 700 000 N/itk-', even at modest applied 
loads. These high pressures result in elastic deformation of materials such 
that the elliptical contact areas are formed for oil film generation and 
load support. The significance of the high contact pressures is that they 
result in a considerable increase in fluid viscosity. Inasmuch as viscosity 
is a measure of a fluid's resistance to flow, this increase greatly enhances 
the lubricant's ability to support load without being squeezed out of the 
contact zone. 

Figure 1 shows the nonconformal surfaces of a ball bearing. The ball 
and race conform to some degree in the section shown in figure 1(a), but the 
section view shown in figure 1(b) clearly exhibits little conformity. We 
are concerned only with nonconformal contacts. 


CURVATURE SUM AND DIFFERENCE 

The undeformed geometry of contacting solids can be represented in gen- 
eral terms by two ellipsoids. The two solids with different radii of curva- 
ture in a pair of principal planes (x and y) passing through the contact 
between the solids make contact at a single point under the condition of 
zero applied load. Such a condition is called point contact and is shown in 
figure 2, where the radii of curvature are denoted by r's. It is assumed 
throughout the lecture that cor v a surfaces, as shown in figure 2, exhibit 
positive curvature and concave < jrf aces, negative curva*'<re. Therefore, if 
the center of curvature lies within the solid, the radiu- of curvature is 
positive; if the center of curvature lies outside the solid, the radius of 
curvature is negative. It is important to note that if coordinates x and 
y are chosen such that 


JL + _L > _L + _L 

r ax r bx r ay r by 


(1) 

coordinate x then determines the direction of the semiminor axis of the 
contact area when a load is applied and y, the direction of the semimajor 
axis. 

The example of a ball bearing is used to illustrate the definition of 
the radii of curvature. A cross section of a thrust-loaded ball bearing 
operating at a contact angle B is shown in figure 3. Equivalent radii of 
curvature for both inner- and outer-race contacts in, and normal to, the 
direction of rolling can be calculated from this figure. The radii of cur- 
vature for the ball - inner-race contact arc 
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The radii of curvature for the ball - outer-race contact 

are 
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Note that the ball - inner-race and ball - outer-race contact inequality (1) 
is satisfied and that the sign convention mentioned earlier has been adopted. 

The curvature sum and difference, which are quantities of some impor- 
tance in the analysis of contact stresses and deformations, are 
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Equations (10) and (11) effectively redefine the problem of two ellipsoidal 
solids approaching one another in terms of an equivalent ellipsoidal solid 
of radii R x and R^ approaching a plane in the manner shown in fig- 
ure 4. From the radius-of-curvature expressions, the radii R x and 
Ry for the contact example discussed earlier can be written for the 
ball - inner-race contact as 
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and for the ball - outer-race contact as 
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GEOMETRICAL separation of ELLIPSOIDAL solids 

The geometrical separation (S ax + Sbx) or (S^y + Sby) between 
two ellipsoidal solids is thus made equivalent to that between a single 
ellipsoidal solid and a plane in the manner shown in figure 4. The geomet- 
rical requirement is simply that at any value of x and y (fig. 4(a)) the 
geometrical separation between the two ellipsoids must be the same as the 
separation between the equivalent ellipsoid and a plane at the same values 
of x and y (fig. 4(b)). From figure 4(a) the following can be written: 


(16) 
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But 2r ax » S ax , so (17) becomes 




(18) 


This is the well-known parabolic approximation to the circular section of 
the solid and is valid as long as the separation is much smaller than the 
radius of curvature. Similar expressions for S a y, S^, and Shy can 
be written, and the expression for the total separation of an ellipsoidal 
solid and a plane (fig. 4(b)) can thus be written as 
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SURFACE STRESSES AND DEFORMATIONS 

When an elastic solid is subjected to a load, stresses are produced 
that increase as the load is increased. These stresses are associated with 
deformations, which are defined by strains. Unique relationships exist be- 
tween stresses and their corresponding strains. For elastic solids the 
stresses are linearly related to the strains, with the constant of propor- 
tionality being an elastic constant that adopts different values for differ- 
ent materials. Thus a simple tensile load applied to a bar produces a 
stress oj and a strain c^, where 
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and 


°1 

E « — * Elastic constant or modulus of elasticity (22) 

e l 

Although no stress acts transversely to the axial direction, there will 
nevertheless be dimensional changes in that direction such that as a bar 
extends axially, it contracts transversely. The transverse strains e 2 
are related to the axial strains ej by Poisson's ratio v such that* 
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where the negative sign means that the transverse strain will be of the 
opposite sign to the axial strain. The modulus of elasticity and Poisson's 
ratio are two important parameters used to describe the material In the 
analysis of contacting solids. 

As the stresses increase within the material, elastic behavior is re- 
placed bj plastic flow, in which the material is permanently deformed. The 
stress state at which the transition from elastic to plastic behavior oc- 
curs, known as the yield stress, has a definite value for a given material 
at a given temperature. In these lectures elastic behavior alone is con- 
sidered. 

When two elastic solids are brought together under a load, a contact 
area develops, the shape and size of which depend on the applied load, the 
elastic properties of the materials, and the curvatures of the surfaces. 

When the two solids shown in figure 2 have a normal load applied to them, 
the shape of the contact area is elliptical, with a being the semimajor 
and b the semiminor axis. It has been common to refer to elliptical con- 
tacts as point contacts, but since these lectures deal mainly with loaded 
contacts, the term el liptical contact is adopted. For the special case 
where r ax « r a y and rt> x ■ r b y, the resulting contact is a circle 
rather than an ellipse. Where r a and r b y are both infinite, the 
initial line contact develops into a rectangle when load is applied. 

Figure 5 shows the contact ellipses obtained with either a radial or a 
thrust load for the ball - inner-race and ball - outer-race contacts in a 
ball bearing. These lectures are concerned with the conjunctions between 
solids with contact areas ranging from circular to rectangular. Inasmuch as 
the size and shape of these contact areas are highly significant to the suc- 
cessful operation of machine elements, it is important to understand their 
characteristics. 

Hertz (1881) considered the stresses and deformations in two perfectly 
smooth, ellipsoidal, contacting elastic solids much like those shown in fig- 
ure 2. His application of the classical theory of elasticity to this prob- 
lem forms the basis of stress calculation for machine elements such as ball 
and roller bearings, gears, seals, and cams. The following assumptions were 
made by Hertz (1881): 

(1) The materials are homogeneous and the yield stress is not exceeded. 

(2) No tangential forces are induced between the solids. 

(3) Contact is limited to a small portion of the surface, such that the 
dimensions of the contact region are small compared with the radii of the 
ellipsoids. 

(4) The solids are at rest and in equilibrium (steady state). 

Making use of these assumptions, Hertz (1881) was able to obtain the 

following expression for the pressure within the ellipsoidal contact shown 
in figure 6. 
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If the pressure is integrated over the contact area, it is found that 


(24) 


p max “ 77Z6 

Equation (24) determines the distribution of pressure or compressive stress 
on the common interface; it is clearly a maximum at the center of the con- 
tact and decreases to zero at the periphery. 
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The ellipticity parameter k can be written in terms of the semimajor 
and semiminor axes of the contact ellipse as 


a 

E 


(* 6 ) 


Harris (1966) has shown that the ellipticity parameter can be written as a 
transcendental equation relating the curvature difference (9) and the ellip- 
tic Integrals of the first S' and second / kinds as 
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A one-point iteration method that was adopted by Hamrock and Anderson (1973) 
car, oe used to obtain the ellipticity parameter, where 

k n-l * k n (30 > 

The iteration process is normally continued until k n+ i differs from 
k p by less than lxlO -7 . Note that the ellipticity parameter is a func- 
tion of the radii of curvature of the solids only: 

k - f(r ,r. ,r ,r. ) (31) 

' ax* bx* ay* by 7 ' 

That is, as the load increases, the semimajor and semiminor axes of the con- 
tact ellipse increase proportionately to each other, so the ellipticity 
parameter remains constant. 

Figure 7 shows the ellipticity parameter and the elliptic integrals of 
the first and second kinds for a range of the curvature ratio 
usually encountered in concentrated contacts. 

When the ellipicity parameter k, the normal applied load F, Poisson's 
ratio v, and the modulus of elasticity E of the contacting solids are 
known, the semimajor and semiminor axes of the contact ellipse and the maxi- 
mum deformation at the center of the contact can be written from the analy- 
sis of Hertz (1881) as 
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where 
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In these equations a and b are proportional to F^-/ 3 an( j 6 p ro _ 
portional to f2/ 3_ 


SUBSURFACE STRESSES 

Fatigue cracks usually start at a certain depth below the surface in 
planes parallel to the direction of rolling. Because of this, special 
attention must be given to the shear stress amplitude occurring in this 
plane. Furthermore a maximum snear stress is reached at a certain depth 
below the surface. The analysis used by Lundberg and Palmgren (1947) will 
be used to define this stress. 

The stresses are referred to a rectangular coordinate system with its 
origin at the center of the contact, its z-axis coinciding with the interior 
normal of the body considered, its x-axis in the direction of rolling, and 
its y-axis in the direction perpendicular to the rolling direction, in the 
analysis that follows it is assumed that y * 0. 

From Lundberg and Palmgren (1947) the following equations can be 
written: 


3F cos t> sin P sin y 
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stress amplitude is defined as 
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For the point of maximum shear stress 
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The position of the maximum point is determined by 
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Furthermore the magnitude of the maximum shear stress is given by 
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SIMPLIFIED SOLUTION FOR STRESSES AND DEFORMATIONS 

The classical Hertzian solution presented in the previous section re- 
quires the calculation of the eilipticity parameter k and the complete 
elliptic integrals of the first and second kinds jr and /. This entails 
finding a solution to a transcendental equation relating k, and / to 

the geometry of the contacting solids, as expressed in (27). This is usu- 
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ally accomplished by some iterative numerical procedure, as described by 
Hamrock and Anderson (1973), or with the aid of charts, as shown by Jones 
(1946). 

Brewe and Hamrock (1977) used a linear regression by the method of 
least squares to obtain simplified equations for k, and /. That is, 
for given sets of pairs of data, {L kj , (Ry/R x )j], J « 1 , 2, .... n}, 
a power fit using a linear regression by the method of least squares re- 
sulted in the following equation: 


IT 


1.0339 



(47) 


The asymptotic behavior of £ and f was suggestive of the functional 
dependence that T and W might exhibit. As a result a logarithmic and 
an inverse curve fit were tried for £ and respectively. The follow- 
ing expressions from Brewe and Hamrock (1977) provide an excellent approxi- 
mation to the relationships between /, JT and Ry/R x : 


S- 1.0003 


0.5968 

or 

y x 


(48) 


jr, 1.5277 ♦ 0.6023 In (R /R ) (49) 

Jr A 

Values of k, £, and W are presented in table 1 and compared with the 
Hamrock and Anderson (1973) numerically determined values of k, /, and 
The agreement is good. _ _ _ 

Using these simplified expressions for k, /, and f and equation 
(34) gives the deformation at the center of the contact 
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(50) 


where 


K - *TTE' ( ■ R Lj \ 12 (51) 

\4.5 r 6 ) 

Note that the load-deflection constant K is a function of the geometry And 
the material properties. 

The results of comparing a with a are also shown in table 1. The 
agreement is again quite good. Therefore the deformation at the center of 
the contact can be obtained directly from (47) to (51). This valuable ap- 
proximation eliminates the need to use curve fitting, charts, or numerical 
methods. 

Figure 8 shows three different degrees of ball-contact conformity: 
a ball on a ball, a ball on a plane, and a ball - outer-ring contact. 

Table 2 uses this figure to show how the degree of conformity affects the 

contact parameters. The table shows that k is not exactly equal to unity 

for the ball-on-ball and ball-on-plane situations because of the approxima- 

tion represented by equation (47). The diameter of the balls is the same 
tnroughout, and the material of the solids is steel. The ball - outer-ring 
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contact Is representative of a 20 y radial ball bearing. A 4.45-N (1-lbf) 
normal load has been considered for each situation. The maximum pressure 
decreases significantly as the curvature of the mating surface approaches 
that of the ball. Table 2 shows that the curvature of the mating surfaces 
Is very Important In relation to the magnitude of the maximum pressure or 
surface stress produced. A ball and ring of high conformity are thus 
desirable from the standpoint of minimizing the stress. 

Table 2 also shows that the area of the contact «ab increases with 
the conformity of the contacting solids. Although this effect minimizes 
contact stresses. It can have an undesirable effect on the force of friction 
since friction force Increases as the contact area, and hence the area of 
the sheared lubricant. Increase in a bearing operating under elastohydro- 
dynamic conditions. The curvatures of the bearing races are therefore gen- 
erally a compromise that takes Into consideration the stress, load capacity, 
and friction characteristics of the bearing. 

In equations (42) to (46) the location and magnitude of the maximum sub- 
surface shear stress are written as functions of t a , an auxiliary parame- 
ter. Furthermore in equation (41) the elliptlclty parameter is written as a 
function of t 4 . The range for 1/k Is 0 £ 1/k < 1, and the correspond- 
ing range for t a is 0 < t d < (1 ♦ vT7)/47 A lTnear regression by the 
method of least squares was used to obtain a simplified formula for t A 
in terms of k, the elliptlcity parameter. That is, for given sets of pairs 
of data {[(1 - t a )j, (1/k), j, j - 1, 2, ..., n], a power fit using a 
linear regression By the method of least squares resulted in the following 
equation: 


I 


a 


0.3044 



l.«55y 


(52) 


The agreement between tnis approximate equation and the exact solution Is 
within *2 percent. Using equation (52) greatly simplifies finding the val- 
ues for the location and magnitude of the maximum subsurface shear stress 
expressed in equations (42) to (46). 


CONCLUDING REMARKS 

An assorted group of topics has been presented dealing with conformal 
and nonconformal surfaces, curvature sum and difference, and surface and 
subsurface stresses in elliptical contacts. Load-def lection relationships 
have been developed for any type of contact. The deformation within the 
contact is a function, among other things, of the ellipticity parameter and 
elliptic Integrals of the first ana second kinds. Simplified expressions 
have been written that allow deformation to be calculated quickly and, gen- 
erally, with adequate accuracy. A simplified expression was also developed 
that enables the location and magnitude of the maximum subsurface shear 
stress to be written direct 1 y without relying on an iterative procedure. 
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APPENDIX - SYMBOLS 


a 

b 

d 

d e 

E 


£' 


£ 

f 

r 

K 

k 

P 

Pmax 

R 

r 

5 

l a 

*.y.z 

e 

r 

6 

‘i 

V 

°1 

T 

T o 


semimajor axis of contact ellipse, m 
semiminor axis of contact ellipse, m 
ball diameter, m 


pitch diameter, m 
modulus of elasticity, N/m* 

effective modulus of elasticity, 




elliptic Integral of second kind with modulus (1 - i/k*) A ^ 
normal applied load, N 

elliptic integral of first kind with modulus (i - i/k*)^* 
load-def lection constant defined in eq. (51) 
ellipticity parameter, a/P 
pressure, N/m* 

maximum pressure within contact, 3 F/*«ab, N/m* 

curvature sum, m 

race curvature radius, m 

geometrical separation of ellipsoidal solids, m 

auxi 1 iary parameter 

coordinate system 

contact angle 

curvature difference 

elastic deformation, m 

strain in axial direction 

strain in transverse direction 

Poisson's ratio 

stress in axial direction, N/m* 
shear stress, N/m^ 


maximum subsurface stress, N/m^ 


Subscripts: 


a solid a 

b solid b 

i inner race 

o outer race 

x.y.z coordinate system 

Superscript: 


( ) approximate 
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TABLE 1. - COMPARISON OF APPROXIMATE AND EXACT FORMULAS 


Radius 


ElUpUcity 


Elliptic Integrals 


Deformation at Center of Contact 

Hallo, 

y\ 

k 

k 

Per- 

cent 

Error 


Second 

kind 



kind 


Numerical, 

From 

Per- 




S 

7 Per- 
cent 
Error 

r 

i 

Per- 

cent 

Error 


Equation, 

6 

Error 

1.000 

1.00 

1.03 

3.00 

1.57 

1.60 1 

.91 

1.57 

1.53 

-2.55 

1.230*10'* 

l.l«8*10' 4 

-5.04 

2.820 

1.99 

2.00 

.50 

1.21 

1.21 C 


2.15 

2.15 

0 

1.020 

1.017 

-.29 

5.314 

3.01 

3.00 

-.33 

1.11 

1.11 


2.53 

2.53 

0 

.897 

.699 

.22 

8.330 

4.01 

4.00 

-.25 

1.07 

1.07 


2.80 

2.60 

0 

.614 

.816 

.25 

11.805 

4.99 

5.00 

.20 

1.05 

1.05 


3.02 

3.01 

-.29 

.756 

.752 

-.53 

15.697 

5.97 

6.00 

.50 

1.04 

1.04 


3.19 

3.16 

-.25 

.706 

.701 

-.71 

19.971 

6,92 

7.00 

1.16 

1.03 

1.03 


3.33 

3.33 

0 

.667 

.662 

-.75 

24.605 

7.87 

6.00 

1.65 

1.02 

1.02 


3.46 

3.45 

-.24 

.636 

.628 

-1.26 

29.576 

8. B0 

9.00 

2.27 

1.02 

1.02 


3.57 

3.56 

-.22 

.608 

.598 

-1.64 

34.869 

9.72 

10.00 

2.88 

1.02 

1.02 1 


3.67 

3.66 

-.25 

.584 

.571 

-2.23 


TABLE 2. - EFFECT OF DECREE OF CONFORMITY 
CONTACT PARAMETERS 

(}For *11 three cases, E’ ■ 2.28X10 11 N/m 2 (steel on steel), 

F - 4.45 N, end r u = r #y * 6.35 mm. For the bell - 
outer ring conUct, d. * 65 mm, 0 » 0°, L » 0.52 
(assume 209 redial ball bearing).] 


Contact 

Parameter* 

Ball on Ball 

IP.ll on Plane 

Ball - Outer 
Ring Contact 

r b« 

6. 35 mm 

- 

-38.9 mm 

r by 

6.35 mm 

* 

-6.60 mm 

R 

1.59 mm 

3.18 mm 

7.26 mm 

yy 

1 

1 

22.1 

k 

1.02 

1.03 

7.330 

7 

1.60 

1.60 

1.03 

j 

1.53 

1.53 

3.38 

a 

0.0465 mm 

0.0586 mm 

0.247 mm 

b 

0.0451 mm 

0.0569 mm 

0.0336 mm 

sab 

0.659*10" 2 mm 2 

1.0«*10' 2 mm 2 

2. 60*10" 2 mm 2 

b 

6. 31*10” 4 mm 

4.87*10“ 4 mm 

2. 56*10 ^ r m 

Pmix 

1.01*10 9 N/m 2 

0.857*10 9 N/m 2 

0.25«*10* N/m 2 










o 4 8 12 16 20 24 28 32 

Curvature ratio, Ry/R, 


Figure 7. - Elliptldty parameter and elliptic Integrals of first and second kinds as 
a function of curvature ratio. 



(c) Bail - outer-ring contact. 
Figure 8. - Three degree? nf cemiurmity. 



